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The 2009 outburst of accreting millisecond pulsar IGR J17511-3057 
as observed by SWIFT and RXTE 
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ABSTRACT 

The twelfth accretion-powered millisecond pulsar, IGR J175 11-3057, was discovered in 
September 2009. In this work we study its spectral and timing properties during the 2009 
outburst based on Swift and RXTE data. Our spectral analysis of the source indicates only 
slight spectral shape evolution during the entire outburst. The equivalent width of the iron line 
and the apparent area of the blackbody emission associated with the hotspot at the stellar sur- 
face both decrease significantly during the outburst. This is consistent with a gradual receding 
of the accretion disc as the accretion rate drops. The pulse profile analysis shows absence of 
dramatic shape evolution with a moderate decrease in pulse amplitude. This behaviour might 
result from a movement of the accretion column footprint towards the magnetic pole as the 
disc retreats. The time lag between the soft and the hard energy pulses increase by a factor 
of two during the outburst. A physical displacement of the centroid of the accretion shock 
relative to the blackbody spot or changes in the emissivity pattern of the Comptonization 
component related to the variations of the accretion column structure could cause this evo- 
lution. We have found that IGR J17511-3057 demonstrates outburst stages similar to those 
seen in SAX J1808.4-3658. A transition from the "slow decay" into the "rapid drop" stage, 
associated with the dramatic flux decrease, is also accompanied by a pulse phase shift which 
could result from an appearance of the secondary spot due to the increasing inner disc radius. 
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1 INTRODUCTION 

Accretion-powered millisecond pulsars (AMPs) are neutron stars 
(NS) that experience transient accretion episodes and show mil- 
lisecond pulsations corresponding to the stellar rotational period. 
Presently, 13 objects of this cl a ss are known (see reviews by 
IPoutanenI I2OO6I : Iwiinandj [20o3 : IPatrund I20T0I) . The pulsations 
arise because the NS magnetic field channels the accretion flow on 
to the NS magnetic poles. Such an accretion flow produces close to 
sinusoidal pulse profiles in most AMPs, but for several of them a 
strong evol ution and peculiar double-peaked pulse shapes are ob- 
served (e.g. lHartman et al.ll2008tllbragimov & Poutanenll2009h . 

The energy spectra of AMPs contain 0.5-1 keV blackbody 
emission from the neutron star surface and a Comptonization 
component dominating at highe r energies and associated proba- 
bly with the accretion shock ( e.g., ^Gierlin ski. Done. & Barre li2002l : 
IPoutanen & Gierlinskillioosi ; [palanga et al.l l2005bl. l2007h . In addi- 
tion, emission from the accretion disc at soft energies (< 2 keV) has 
also been detected in XMM-Newton observations of XTE J1751- 
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305 jGieriiirski & Poutanen |l2005^■ XTE JI807-294 jpalanga et all 
l2005al) , SAX J1808.4-3658 (Patruno et al. lOoj) and receii tly in 
IGR J17511-3057 (hereafter IGR17511, IPapitto et aljboiol) . The 
reflecti on component of moderate amplitude is also detected i n 
AMPs jGieriiilski & Poutanenl200^ ; llbragimov & Poutanerj2009h . 

The pulse profi les are often rather sinusoidal with a slight 
skewness jPoutanerJ [ 2OO 6I and show clear energy dependence. 
Pulses at soft energies peak at a later ph ase resulting in "soft lags", 
as firs t was seen in SAX J 1808. 4-365 8 jCui. Morgan. & TitarchukI 
Il998l) . The origin of the soft lag is most likely related to the 
different angular emission patterns of the blackbody and Comp- 
tonized components, which naturally explains why the phase lags 
seem to saturate at ~ 7-10 keV, where the emission of the 
blackbody component becomes negligible llGierliiiski et alj|2002l : 
IPoutanen & Gieriiirskillioosi : iGieriiiiski & Poutanenll2005h . Some 
AMPs however show also a (not fully understood) decrease in the 
lag above ~ 10 keV (see n in IGR 10029 1-1-59 34 as indicated by 
iGallowav et al] |2005| and iFalanga et all l2005bl: IGR175 11 migh t 
also have this decrease above ~ 20 keV, see Falanga et ani201 1[) . 
Information about this lag can be used to study the properties of the 
accretion shock and the structure of the hotspot. 



2 A. Ibragimov, J. J. E. Kajava and J. Poutanen 



1.1 IGR J17511-3057 

IGR17511 was discovered on 2009 September 12 (MJD 55087) by 
INTEGRAL observatory jBaldovin et"aLll2009|) duri ng the Galac- 
tic bulge monitoring program dKuulkers et al.'7007^. The 245 Hz 
pulsations were detected by RXTE ( Markwardt et al. 20093) con- 
firming the AMP nature of IGR17511. A Chandm/HETG obser- 
vation provided the source position of (J2000) RA=17''51 '"08':66, 
Dec=-30°57'41'.'0 (la error of 0"6. iNowak etal]|2009l) . A near 
infrared counterpart of magn itude Ks = 18™0 ± 0.1 was iden- 
tified by [Torres et"al1 (|2009|) within the Chandra error box, but 
no r adio counterpart was detected with a 3 cr upper limit of 0.10 
mjy jMUler- Jones. Russell. & Migliaril2009l) . The source faded be- 
yond RXTE detection l imit after 2009 October 11 (MJD 55113, 
iMarkwardt et al.ll2009ah . 

Type I X-ray bursts were observed in IGR 175 11 with Swift 
( IBozzo et alj 20091) and bur st oscillations imme diately after with 
RXTE l lWatts etal.l2009l see lFalanga etai]|201ll for the analysis of 
all detected bursts). Several distance constraints have b een reported 
based on these data. The analysis of Swift data by IBozzo et alj 
( I2OIOI) yielded an upper limit on the distance of 10.1 ± 0.5 kpc, 
derived using the empirical relation of the Eddington limit Lodd ~ 
(3.79 ± 0.1 5) X 10'^^ erg s~^ fo r the photospheric radius expan- 
sion bursts dKuulkers et al.ll2003h. Another upper limit of 7.5 kpc 
was derived by Falanga et alj ^ 20 111) v ia the independent a nalysis of 
type I bursts. Using the s ame method, iPapitto et al. reported 
a similar upper limit as IBozzo et al 



Table 1. Data groupings 



Group code MJD interval/RXTE MJD interval/SW/FT 



( 1201 



from XMM-Newton 



data, but the analysis of RXTE data by Altamirano et alj _([201(lf) 
gave a tighter upper limit of 6.9 kpc. Altamirano et al. ( 201(1) also 
used the distance approximation of iGallowav et al.l ( |2008|) . that re- 
sulted in an upper limit of 4.4 kpc for a NS of mass 1.4Mq , radius 
10 km and hydrogen mass fraction X — 0.7. The corresponding 
upper limit for X = would be 5.76 kpc. but absence of hydrogen 
is inconsistent with the fact that the companion of IGR 175 11 seem s 
to be a main sequence star dPapittoetalfcOld : iRiggio et al.ll201 1[) . 
In light of these constraints, we adopt the distance value of 5 kpc. 

The source light curve shows an exponential flux decay, 
commonly seen in A MPs (a "slow decay" in the terminology of 
iHartman et al ]|2008h . The pulse profiles are single-peaked, indicat- 
ing that we probably see the emission coming mainly from one 
emitting spot on the neutron star surface (i.e. contribution from the 
secondary spot does not produce a distinct secondary feature; how- 
ever, "flattened" pulse profile minima may suggest a presence of 
secondary spot emission). 

In this paper, we present the results of our spectral and timing 
analysis of IGR17511 based on Swift and RXTE observations. We 
study the evolution of phase averaged- and phase resolved spectra, 
phase lags and pulse profile changes during the outburst. 



2 OBSERVATIONAL DATA 

The RXTE data covering the outburst of IGR 1751 1 (ObsID 94041) 
were reduced using HEASOFT 6.8 and the CALDB. We used data 
taken both by RXTE/PCA (3-25 keV) and HEXTE (25-200 keV). 
Standard 0.5 per ce nt systematic was applied to the PCA spectra 
djahoda et alj|200d) . To keep the calibration uniform, we used data 
from PCA unit 2 only. T he source spectrum is c ontaminated by the 
Galactic ridge emission dRevnivtsev et alj2009h . To take it into ac- 
count for RXTEfPCA spectra, we have produced a spectrum from 
the observations, where both IGR17511 and nearby AMP XTE 
J1751-305 were in a quiescent state (MJD 55115 - 55126). This 
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55087.9- 


-55109.3 


55087.8-55107.5 
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55087.9- 


-55088.9 


55087.8-55088.7 
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55089.2- 


-55090.5 


55089.6-55090.8 
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55091.2- 


-55094.0 


55092.8-55093.9 
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55094.0- 


-55096.8 


55094.4-55095.5 
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55097.2- 


-55099.5 




6 


55100.3- 


-55101.9 




7 


55102.2- 


-55104.9 


55102.4-55105.0 


8 


55105.4- 


-55109.3 


55107.0-55107.5 



spectrum, mainly affecting channels below 15 keV, was subtracted 
from all spectral and timing data. 

For Swift/XKT observations, we only considered window- 
timing mode data, because the photon-countin g mode data suf- 
fered from photon pile-up dBozzo et al. I I2OIOI) . We reduced the 
Swift/XKT data with XRTPIPELINE v.0.12.3 using standard filter- 
ing and screening criteria for the event selection. We used circular 
regions of 20 pixel radius to extract the spectral data. The XRT- 
EXPOMAP task was used to generate the exposure maps and the 
ancillary response files were generated with the XRTMKARF task to 
account for different extraction regions, vignetting and point spread 
function corrections. Ancillary response files (ARFs) of individ- 
ual Swift/XKT snapshots were averaged together; each ARE was 
co-added with "weight" equal to the relative contribution of pho- 
tons detected in the snapsot to the overall photon number collected 
from all snapshots. The Swift/XKT redistribution matrices (v.Oll) 
were taken from the CALDB. Swift/XKT and RXTE/YIEXJE spec- 
tra were grouped such as eac h bin contained at least 200 counts. 

The type I X-ray bursts dBozzo et al.l2010l;|Papitto et al.[2010l ; 
lAltamirano etai] I2OI0I ; iRiggio et alj 1201 ll ; lEalanga et alj hoTlh 
were screened out from our analysis. T he spectral analysis was 
done using XSPEC V.12 dAmaudll 19961) . Uncertainties of spectral 



and timing best-fitting parameters correspond to the 90 per cent 
confidence level, unless otherwise stated. 



3 SPECTRAL ANALYSIS 

In this section, we describe the results of our spectral analysis of the 
source. In order to improve statistics, we have grouped individual 
spectra as described in Table [T] 



3.1 Spectral model 

The spectrum of IGR1751 1 is typical for an AMP and can be de- 
scribed as a composition of accretion disc emission (around 1-2 
keV, not visible in the RXTE range), blackbody originating from 
the hotspot (2-10 keV) and hard X-ray tail generated by thermal 
Comptonization in accretion shock located above the neutron star 
surface (contributing in the whole range of 1-200 keV and dom- 
inating above 1 keV), similar to other objects of this kind (e.g. 
'Gieriinski et al.'l 2002l; [cieriinski & PoutanenI l2005t lEalanga et alj 
i2005a b, 2007; see Eig. [T]l. To model th ermal Comptonization 
contin uum we used the COMPPS model of IPoutanen & SvenssorJ 
dl996h . A fluorescent iron line at 6.4 keV and the Compton re- 
flection of the COMPPS component dMagdziarz & Zdziarski|[l995l) 
were also included in the fitting. The spectral model includes 
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Figure 1. The joint observed spectrum of IGR175 1 1 from RXTE and Swift 
satellites collected during the entire outburst (group T). Red, green and 
blue data points represent Swift/XRT, RXTE/PCA and RXTE/HEXTE re- 
spectively. Solid, dotted, long-dashed, dashed and dot-dashed curves repre- 
sent the total spectrum, Comptonization continuum, reflection and iron line, 
blackbody emission from the hotspot and the disc blackbody, respectively. 
Lower panel show the residuals of the fit. The fit parameters are given in 
Table |2] (group T for joint RXTE and Swift spectra). En'or bars con'espond 

to 1(T. 



also interstellar absorption model PHABS parametrized by the hy- 
drogen column density A^'h. The described approach corresponds 
to PHABS*(DISKBB-FBBODYRAD-F COMPPS-I-DISKLINE) model in 
XSPEC. 

For COMPPS we assumed a slab geometry and the model is 
characterized by the Thomson optical depth tt and the temperature 
of the hot electrons To. The seed photons for Comptonization have 
a temperature Tscod and the surface area is denoted as Eshock- The 
blackbody component has a temperature Tbb and its surface area is 
denoted as Espot- The apparent spot radii at infinity can be com- 
puted from the BBODYRAD and COMPPS model normalizations: 
E — nB? — iiKD\q, where K is the normalization obtained 
from fits, R is the apparent radius in kilometres and Dio is the dis- 
tance in units of 10 kpc. We denote these radii as i?spot and iishock 
for the BBODYRAD and COMPPS, respectively. The Compton re- 
flection from the accretion disc is parametrized by the amplitude 
K = Q./ 2ti, where n is the solid ang le covered by the reflecting 
medium jMagdziarz & Zdziarskil 19951) . We used DISKLINE model 
to model the iron line and we fixed the inner disc radius rin = lOrg, 
where r-g = IGMji? is the Schwarzschild radius. We also as- 
sumed that the radial emissivity profile for the illuminating contin- 
uum flux cx r"^. For the DISKLINE, as wefl as for COMPPS, we 
assumed an inclination of i = 60°. The DISKBB model component 
was only included in the joint analysis of Swift and RXTE data (Sec- 
tion l3.2t to account for the soft X-ray accretion disc emission. The 
respective model parameters are the inner disc temperature Tin and 
the apparent inner disc radius J?in, which can be computed from 
the model normalization Kdhh as i?in = D\o \J -ft^dbb / cos i. 



3.2 Phase-averaged spectra from RXTE and Swift 

We used the SwiftKKl data (in range 0.6-8.0 keV) to study the 
soft X-ray emission of IGR 17511. Initially, we fitted group T to 
constrain the disc parameters and the absorption column. The emit- 
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Figure 2. The best-fitting parameters for the Comptonization-based model 
of Section [331 The absorption corrected flux in 3-20 keV band is in units 
of 10^^ erg cm^'^ s~^. 



ting areas of the hotspot and the shock were assumed equal (while 
the spectrum is statistically good to fit these areas independently, 
a large number of free parameters would lead to large uncertain- 
ties; in fact, these areas must be similar on physical grounds, see 
llbragimov & PoutanenI |2009|) . The best-fitting results are shown 
in Table [2] in particular we find that the emission below < 3 
keV has cle ar signatures of the accretion disc (see Fig. [T] and also 
IPapitto eTal, 2 010, fig. 9). Fits with free interstellar absorption lead 
to A^H = (O.SSloii) X 10^^cm"^ Tin = 0.24 ± 0.07 keV and 



i?in = 40 



-1-46 



km. 



The aforementioned best fitting values are subjected to sev- 
eral uncertainties. The values of the disc parameters Tin and i?in 
are tightly correlated with the absorption value A^h- In addition, the 
derived value of i?in depends on the assumed distance and it should 
be corrected for two effect s in ord er to obtain a realistic radiu s. As 
discussed in [KubotaeTal] i 19981) and lGierMski et alj l ll999h . the 
derived radius should be larger by the square of the colour correc- 
tion factor /c. This factor h = 1.7 was computed fo r accretion 
discs around black holes bv lShimura & Takaharal ( ll995h . However, 
in the case of AMPs, there is a stark difference in the sense that 
the disc is irradiated by the emission from the hotspot, which casts 
an uncertainty to this value. Furthermore, Rin is also affected by 
a corr ection due to the inner boundary condition dOierliiiski et alj 
1 19991) . but this (of the order of unity) factor is not accurately known 
in the case of accretion onto a magnetized star. Therefore, the value 
of the inner disc radius should be taken as an order of magnitude 
estimate. 

Other sp ectral parameters (se e Table |2) are consistent with 
the findings of lFalanga et alj ( 12011 ). However, there ar e small dif- 
ferences between our results and Papitto et al.l boiol) especially 
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Table 2. Results of spectral fitting with the Comptonization model with constant areas for the blackbody and Comptonization compo- 
nents (Sect. [X2]|3.3) . Left column indicates the data group. Letter "F" indicates a fixed parameter 
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RXTE and Swift, interstellar absorption yield A^h = {O.88IQ-24) x lO^^cm-^. 
i?XT£only. 



Table 3. Results of spectral fitting with the Comptonization model using independent areas of the hotspot blackbody and Comptonized component using RXTE 
data only. Sect. 13. 31 Left column indicates the data group. 



Group Tc TT T.eed K Tbb EW Espot S.hock X^/d.o.f 

keV keV keV eV km^ km^ 



T 33I3 1.61 ±0.18 l.lSloa? 0.45j;oi2 0.60 ±0.13 46l|* 63l26* 34^2? 190/160 

1^ 35 ±3 1.54l0;l« 1.161°;?^ 0A9tt\t 0.6lt°\l 6ll|? SOtH'' 431^8 148/I6O 
5-8 3ltl l-66lE!;29 l-SltoJg OAGtoll 0.63to\l < HO 441?? 23^^^ 206/160 



for the values of tt and Tsccd- The most likely reason for these 
differences is because the high energy cutoff cannot be accu- 
rately determined in short HEXTE exposures. Also, the cross- 
cali bration issues between SwiftfX.KT and XMM/EPIC instruments 
(see lTsuiimoto et al]|201lh might cause the differences. 

In further analysis, we found that the Swift/XKF data do not 
have enough statistics to reliably constrain the disc component for 
individual fits of groups 1-8. Therefore, we could not look for 
changes in the disc parameters during the outburst and in the fol- 
lowing sections we only consider RXTE data when we study the 
evolution of the spectral parameters. 



3.3 Phase-averaged spectra from RXTE 

We began our RXTE-only spectral analysis by fitting the model 
with independent emitting areas for the hotspot and the shock (free 
blackbody and Comptonized component normalizations). Because 
the accretion disc does not contribute to the flux in the PCA band 3- 
25 keV, we omit the DISKBB component from the following spec- 
tral fits and use our best fitting value of A'h = 0.88 x 10^^ cm"^ 
in the fitting of iMiTB-only data. The spectrum for group T and 
"joined" groups l-^ and 5-8 have sufficient statistics to fit these 
areas separately. 

The results of the fitting are shown in Table|3] We find that the 
ratio between the areas is compatible with constant, although the 
statistical errors are large to ma ke a firm conc l usion . We note that 
analysis of the same source by IPapittoetal.1 Jioioh and of SAX 



J 1808. 4-365 8 by llbragimov & Poutanei] | |2009|) both suggest, in 
agreement with our result, that the blackbody area is 2-3 times 
larger than the area of Comptonized emission. These fits indicate 
a decrease of the emitting areas, as expected from a gradually in - 
creasing inner disc radius dPoutanen. Ibragimov. & Annalall200^ . 
The simultaneous decrease of the iron line equivalent width sup- 
ports the expected physical picture (note that the reflection ampli- 
tude should decrease as well, but observation uncertainties do not 
allow us to constrain it reliably). 

The group T and the "joined" groups 1-4 and 5-8 also al- 
low for independent fitting of Tc and reflection amplitude 5R. The 
best fitting parameters are shown in Table |2] We note that the ac- 
tual value of 5R is subject to the chosen continuum model (e.g., 
iDone. Gierliiiski. & Kubot3l 20071) . 

Spectra for individual groups 1-8 do not have enough statis- 
tics to constrain Te and (that turns out to be uncertain and com- 
patible with zero). Therefore, we fixed ro=30 keV and — 0.3 
(as found in the Swift and RXTE fits. Sect. 13.21 ) and obtained fits 
for individual data groups. Furthermore, we adopted equal emit- 
ting areas for the Comptonization and blackbody components for 
these groups, because the data quality does not allow us to fit them 
independently and they should be similar on physical grou nds (see 
iGierlinski & Poutanenllioosi : llbragimov & Poutanenll2009l) . 

The time evolution of the spectral parameters is shown in Fig. 
[2] and in Table |2l The optical depth is initially roughly constant 
Tt ~ 1.85 for groups 1-7, but later it drops to tt = 1.74 ± 0.03 
for group 8. This shows (in connection with fixed Tc value) that 
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Figure 3. Results of fitting the per-orbit pulse profiles in (a) 2.1-3.7 keV and (b) 9.7-23.1 keV with expression 0. Top panels: the amplitudes of the 
fundamental ai (squares, red) and overtone 02 (diamonds, blue). Bottom panels: phases of the fundamental (pi (squares, red) and overtone (/>2 (diamonds, 
blue, shifted by phase 0.5 for clarity). Typical error bars are shown for a few points. Note the changes in amplitude and phase of the fundamental around MJD 
55112, addressed in Section lTl] 



the spectrum softens in the end of the outburst. It is also noticeable, 
that Espot, Tsccd and Tbb are decreasing slightly during the out- 
burst. The decrease in Espot is most likely caused by the change 
of Rin in the course of the outburst jPoutanen et al.ll2009l) . When 
the flux drops as the mass accretion rate goes down, iiin increases 
(it is likely proportio nal to the Alfven radius whi ch has a Af^^''^ 
dependence, see e.g. iFrank. King. & Raind[20o3) . Assuming that 
the magnetic field is a dipole, the outer boundary of the hotspot is 
controlled by the current position of Rin and therefore the increase 
in Rin leads to decrease in Espot- 



4 TIMING ANALYSIS 

The pulse shape of an accreting pulsar contains important informa- 
tion abo ut physics of emission and geometrical parameters of the 
system (Poutanen & Gierliiiski 2003, Leahy, Morsink, & Cadead 
2008 ; Poutanen 2008; Ibragimov & Poutanen 2009; Poutanen et al, 
2009h To obtain the pulse profiles, we used the ephemeris of 
Riggio et"^ ( I2OO9I) . In general, the pulse profiles of IGR 175 1 1 are 
single-peaked, close to symmetric and without a prominent sec- 
ondary maxima. 



4.1 Harmonic content evolution 

From our spectral analysis (e.g.. Fig. [T) we see that the hard X-ray 
part of the spectrum above 10 keV is dominated by thermal Comp- 
tonization, that probably takes place in the accretion shock. The soft 
part (below 10 keV) includes blackbody emission from the hotspot. 
Furthermore, around 2 keV there is em ission from the accretion 
disc (see Fig. [Tl and IPapitto et al JI2OI0I fig. 9). Consequently, we 
chose energy ranges of interest as 2.1-3.7 and 9.7-23.1 keV; the 
first interval contains a large fraction of the hotspot blackbody ra- 
diation (and some part of disc emission which is non-pulsating), 
while the second band c ontains only Compto nized emission. While 
a throughout analysis bv lRiggio et al.l ( l20Ilh demonstrated that one 
needs three (and, sometimes, even four) harmonics to fully describe 
wide-energy (2-25 keV) pulses, we find it acceptable to utilize two 
Fourier harmonics for description of pulses in narrow energy bands 
and trace the pulse profile changes. We note that wide-energy pulse 
profile is in fact a superposition of different pulse shapes seen at 
different energies and this might affect the harmonics decomposi- 
tion. However, indeed the comparison between two-harmonics fit 
and the actual pulse profiles reveal some deviations from a smooth 
fit near the pulse profile minima, that can be due to the antipodal 
spot contribution or due to addition al absorption at cer tain phase 
(by e.g. the accretion column), which lRiggio et al1 ( l20I ih modelled 
with a 3rd harmonic (as it follows from their fig. 2). We fitted the 
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MJD Phase Phase 



Figure 4. Left panels: evolution of harmonics content for groups 1-8. Stars and triangles represent 2.1-3.7 keV, diamonds and squares - 9.7-23.1 keV. (a) 
The amplitudes of the fundamental ai (upper points) and overtone 02 (lower points); (b) phases of the fundamental 0i and (c) overtone 4i2 versus time. Right 
panels: pulse profiles for 2.1-3.7 and 9.7-23.1 keV (blue and red histograms, respectively). For amplitudes and phases (panels a-c) the errors are at the 90 
per cent confidence level, while for pulse profiles the errors bars correspond to la. 



pulse profiles collected from each RXTE orbit with the following 
expression 



F{(j)) = F{1 + ai cos[27r((?i - (^1)] + az cos[47r(0 - ^a)]}, 



(1) 



where ai, 02 are amplitudes and (j>\, <j)2 are phases of the funda- 
mental and the overtone, respectively. 

The fitting results for aforementioned energy bands are shown 
in Fig. (3] The amplitude of the fundamental decreases with time 
(and with flux). While the general trend is rather smooth, some ir- 
regularity in fundamental amplitude can be noticed around MJD 
55098, and phase of overtone experiences small shift around 
MJD 55093. After MJD 55100, we observe a "drift" in the 
phase of fundame ntal, seen clearly on Fig. |4] (see also fig. 4 in 
iRiggio et alj|201 fl) . This drift is larger at soft energies, which in 
turn affects the value of the phase lag (see Sect. 14. 2t . The har- 
monic content on soft and hard energies behaves very similarly, 
which is also seen in SAX J1808.4— 3658 during its "slow decay" 
stage jlbragimov & PoutanenI I2OO9I) . After the slow decay, SAX 
J1808.4-3658 has shown quite a different evolution of pulse shapes 
at different energies. For our case, there is a clear evidence that the 
source experienced a likely transition to the "rapid drop" stage (see 
Section |431 ). However, its flux dropped very quickly below the de- 
tection level, making it impossible to study pulse profile evolution 
further. 



4.2 Evolution of phase lags 

In AMPs the pulses in the soft and hard energies do no t arrive in 
phase , but there is an energy-dependent phase lag (e.g.. ICui et alj 
1 19981) . We determined the phase lags by fitting the pulse profiles 
at a given energy with expression ([T} and finding the phase differ- 
ence relative to the reference energy band 2.1-3.7 keV. This way 
we find the phase lags in each harmonic as well as the pulse maxi- 
mum lag coiTesponding to the phase difference between fitted pulse 
maxima. In IGR17511, the phase lag is negative (i.e. pulse peaks 
on soft energies at a later phase) and shows a gradual increase 
from 3 keV to approximately 10 keV, where the lag val ue nearly 
"saturates". This behaviour is typical for AMPs (see e.g. Cui et al 



1998l:lGierlinski et al.ll2002llHartman. Watts. & Chakrabartvll200' 



a 



Falanga et al.l 201 ih . The phase lag of the overtone is poorly con- 
strained: for our groups 1-8 it is noticeable that overtone best- 
fitting value decreases from to ~ —100 /.ts in the energy range 
3-7 keV, and then remains constant or slightly reduces. But in all 
cases it is compatible with zero and is determined with the uncer- 
tainty more than 100/is. 

The most noticeable effect we saw in the data is a gradual in- 
crease of the phase lag (measured between 2.1-3.7 and 15.5-23.1 
keV) from 200 to 400 /is during the outburst, as illustrated on the 
right panel of Fig. |5] Interestingly, fundamental shows a straight- 
forward trend, while pulse maximum lag changes noticeably only 
close to the end of the outburst. A comparison of the pulse profiles 
obtained at various dates indicates that the lag increases because 
the pulse maximum at soft energies (2.1-3.7 keV in our case) shifts 
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Figure 5. Left panels: pulse profile fits with equation {T] for energies 3.3-3.7, 4.5-4.9, 6.5-8.1 and 15.5-23.1 keV (solid, dotted, dashed and dash-dotted lines, 
respectively). Middle panels: pulse maximum lags (stars) and the lags of the fundamental (diamonds). The overtone lags have large eiTors and therefore are not 
shown here. The observational groups are indicated on the panels. Right panel: evolution of the lags between 2.1-3.7 and 15.5-23.1 keV energy bands of the 
pulse maximum and of the fundamental. The pulse maximum lag was computed as the phase difference between the maxima at different energies, determined 
via fitting the observed pulse shape with expression ([T). 
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Figure 6. The pulse profiles (2-60 keV) from MJD 55111.01 (blue trian- 
gles) and MJD 551 12.07 (red asterisks), demonstrating a significant shift of 
the pulse maximum and an abrupt drop in the amplitude (due to coinciden- 
tal outburst of a neighbouring AMP XTE J175 1-305). Blue solid and red 
dashed curves are the respective best-fitting approximations with expression 
(T). Error bars correspond to 1 a. 



to a latter phase, while the maximum of the high-energy pulse (9.7- 
23.1 keV) shifts in parallel, but in a less pronounced way. To illus- 
trate that, in Fig. |4]we plot a set of pulse profiles from a few time 
intervals and the evolution of their harmonic content. 



4.3 Pulse profile changes at MJD 55112 

The pulsations became completely undetectable shortly after MJD 
55112, and the last useful observations do not contain a reliable 
statistics to obtain a well-defined energy-resolved pulse shape. 
However, the shape in the whole RXTEfPCA range (approximately 
2-60 keV) can be used to locate the pulse maximum. In Fig.|6]we 
show two pulse profiles from the adjacent observations centered 
at MJD 55111.01 and MJD 55112.07 together with the respec- 
tive harmonical fits. The relative amplitude abruptly decreases be- 
cause the average X-ray emission is modified by simultaneous out- 
burst of A MP XTE J175 1-305 in the field of view of RXTEfPCA 
([Markward t et al.l2009ah . It is clear that the phase of the pulse max- 
imum has shifted forward by about 0.1-0.2. By analogue with SAX 
J1808.4-3658, we ca n speculate that the source beg an to shift into 
the rapid drop stage jlbragimov & Poutanenll2009l) . As noted by 
iFalanga etal.N201 ij), the exponential trend typical to the slow de- 
cay outburst stage is followed by a faster, linear drop of flux around 
MJD 55107, few days prior to detected pulse evolution. Fig.[3]re- 
veals that in parallel with the linearly decreasing flux trend, the 
phase of fundamental starts to increase slowly, ending up in a sharp 
phase jump at MJD 55112. We can speculate that the accretion disc 
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Table 4. Harmonic fits by expression (T} to the pliase-resolved apparent 
areas. The average apparent areas E (in units of km^) are computed for the 
distance of 5 kpc. 



55090 



55120 



Figure 7. The lightcurve of IGR17511 during the 2009 outburst. The 
flux in the 3-10 keV range is corrected for absorption. Blue stars corre- 
spond to the RXTE/PCA data, and red diamonds are for Swift/XKT (Photon 
Counting mode). Two vertical dotted lines denote observations centered at 
MJD 55511.01 and MJD 55112.07 that are shown in Fig. [6] Although the 
RXTE/PCA observations are contaminated after MJD 55112 by the out- 
burst of a nearby AMP XTE J175 1-305, Swift allows to determine the flux 
of 1GR17511 unambiguously revealing an abrupt flux fall consistent with 
being the start of the "rapid drop" outburst stage. EiTor bars coiTespond to 
la. 



starts to recede from the neutron star and at some point the antipo- 
dal spot appears to our view changing the pulse shape. 

While RXTEfPCA is a non-imaging instrument and it is im- 
possible to separate flux from IGR17511 and XTE J1751-305, 
we were able to estimate flux received from IGR17511 using 
Swift/XKT. The Swift count rate li ghtcurve has been produced using 
online XRT product generatoiQ jEvans et"aL 1 1200^ and converted 
to energy flux in 3-10 keV band using the webPIMMS too0 usmg 
powerlaw with photon index 1.7 that is suitable for our object in 
the mentioned energy interval. We also obtained the RXTE/PCA 
light curve in the same energy range for cross-calibration. The 
resulting lightcurve of the outburst is shown on Fig. [7] A sharp 
abrupt fall of the flux ha s been observed on October 9 (MJD 55113, 
iMarkwardt et al followed by a non-detection of the source 

in the subsequent pointings. This strongly supports the conclusion 
that the source entered the rapid drop stage. Since there is no Swift 
observation around MJD 551 12, it is not possible to reliably correct 
the amplitude of the respective pulse profile shown on Fig.|6] 

4.4 Phase-resolved spectrum 

Phase-resolved spectro scopy was performed f or 1998 outburst of 
SAX J1808.4-3658 bv iGierliiiski et all ilOOA for the 2002 out- 
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burst bv llbragimov & Poutaiienl (l2009h andlWilkinson et all ( |2010|) 
and for XTE J175 1-305 bv iGierlihski & PoutanenI ( |2005[) . In the 
former work it was concluded that the energy dependence of the 
pulse profiles and phase lags can be explained by a simple model 
where only normalizations of the hotspot blackbody and Comp- 
tonization tail vary. Similar to these papers, we have generated the 
phase-resolved spectra for group T and used the phase-averaged 
spectrum (Section |33} as a reference. We fixed all the parameters 
except the blackbody and Comptonization normalizations. The re- 
sults are shown in Fig. [8] The coefficients for expression ([T) that 
describe the modulation in the apparent areas of the blackbody 



Espot and Comptonized tail Eshock are shown in Table |4] Lag be- 
tween components is large: fundamental components have phase 
difference of 0.18, which corresponds to ~730/is, while the lag 
of the overtone is difficult to determine precisely. Phase-resolved 
fits of groups 1-4 and 5-8 suggest slight increase of component 
lags that reflects the changes in the pulse profile. The determined 
effective areas are shifted comparing to Table |3] values, since we 
use continuum shape determined using equal areas of blackbody 
and Comptonized emission (i.e., given in Table |2j. Utilizing con- 
tinuum parameters obtained with independent emitting areas yield 
in area values compatible with those in Table |3] however, the ob- 
served phase difference remains the same. 

It is interesting to note that the measured ~730/is lag between 
components is much larger than the one observed in actual pulsa- 
tions. The effect also seen in similar studies, e.g., of SAX J1808.4- 
3658 (GierUnski et al. 2002; Ibragimov & Poutanen 2009). It is a 
natural consequence of the fact that in phase-resolved spectroscopy 
we decompose the observed joint spectrum into separate physical 
components, while the pulse that we observe in the soft X-ray band 
is a mix from these components (in nearly equal proportion, see 
Fig. [TJ. Modelling of time lags (described below in Section 15.3b 
confirms the fact that the lag between components is indeed much 
larger than the one between the observed pulses at different ener- 
gies. 



5 DISCUSSION 
5.1 Spot size 

The apparent spot size at infinity i?spot can be related to the 
actual (circular) spot radius by taking into account relativistic 
light bending, as descr i bed in j Poutanen & Gierliiiski ( 200^ and 
llbragimov & PoutanenI ([200^|j We can use the emitting areas 
shown in Table |2] for the group T to obtain estimations of the 



^ http://www.swift.ac.uk/user_objects/ 

^ http://heasarc.nasa.gov/Tools/w3pimms.html 



^ Note that the relation is derived for a slowly rotating pulsar and black- 
body emission. 
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Figure 8. Results of the phase-resolved spectral analysis. The best-fitting pai'ameters (except normalizations) are frozen at the values obtained in the corre- 
sponding phase-averaged fit (Sect. l33l Table|2). Left panel: apparent emitting areas of the two components (blue and red curves represent blackbody Espot 
and Comptonization Eghock components, respectively). Right panels: Sspot vs Sghock- The parameters of the fits with a harmonic function are shown in 
Tableg] 



spot size. We take a star with A/*=1.4 M© and R,= 10.3 km 
(2.5 Schwarzschild radii). The actual spot radius depends on the 
system inclination i and the spot colatitude 6 (see sect. 5.1 of 
llbragimov & Poutanenll20()9l) . The smallest and the largest possi- 
ble radii correspond to i — 9 — 0° and i = 6 = 90°, respectively. 
For best-fitting normalizations of the RXTE-only fits, the interval 
of radii is 4.5-7. 1 km, and for the joint RXTE and Swift spectrum it 
is 4.0-6.4 km. 

Besides the uncertainty in the areas because of a weakly con- 
strained distance to the object, the numbers quoted are subjected 
to uncertainty that resides in the unknown relation between the 
emitting area of the blackbody and Comptonization components 
(which we assumed to be equal). For SAX J1808.4-3658, the time- 
average spectrum (for the slow decay outburst stage) suggested that 
the blackbody area is twice as large as the Comptonized one, and 
the likely error in the area is about 50 per cent; additional uncer- 
tainty due to a colour-correction appears if the emission is different 
from a blackbody, however, for atmospheres heated from above this 
correction should not play a significant role JPoutanen & Gierlinskil 
[2003; Ibragimov & Poutanen 2009). Finally, the estimation is valid 
for the filled circular spot, while in reality the spot shape can be 
rather different (see Section l5!2t . 




5 10 15 20 25 30 35 40 
Spot colatitude (degrees) 

Figure 9. Contour plots of the constant oscillation amplitude at the 
plane inclination - spot cola titude. The curves are computed using ex- 
pression (22) in .P outanen & Gierliriskil j2003l) (see also Sect 5.1 in 
llbragimov & Poutanenll2009l) assuming a 1.4Mq neutron star with radius 
of 10.3 km and the spot size of ij^o = 5 km. 



5.2 Oscillation amplitude and geometry 

IPoutanen & Gierlinskil l l2003h have derived an expression that re- 
lates oscillation amplitude of the pulsar to the system inclination 
and spot colatitude for a case of large, blackbody-emitting, always 
visible spot on the surfa ce of a slowly rotating star (s ee equation 
(10) there and sect. 5.2 in llbragimov & Poutanenll2009l) . While our 
hotspot is non-blackbody and the star rotates rapidly (and there are 
arguments against filled circular spot shape, see below), we can 
compare the observed relative amplitude of the fundamental with 
this analytical relation to obtain a zero-order estimate of system's 
geometrical parameters. Taking the apparent spot radius Roo ~ 5 
km, we can obtain a dependence of the amplitude on the inclination 
and spot colatitude which (and typical neutron star parameters) are 
shown on Fig. [9l We remind that in our case the observed ampli- 



tude of the fundamental is about 22 per cent at the beginning of the 
outburst, gradually decreasing to 15 per cent (Fig. |3)- Assuming 
60° inclination we get the spot colatitude of about 15°. A decrease 
in the spot colatitude can cause a corresponding decrease of the 
amplitude. 

What could be a reason for decreasing spot colatitude dur- 
ing the outburst? Present understanding of the AMP geometry is 
that of the inclined dipole with the accretion disc disrupted at some 
truncation radius (likely proportional to the Alfven radius, see e.g. 
iLong et al] |2005|) by the magnetic field. Matter flows along the 
magnetic field lines and falls on the neutron star surface form- 
ing a hotspot. As the accretion rate drops, the Alfven radius in- 
creases and the disc recedes from the neutron star. The matter then 
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accretes along the magnetic field lines that touch the star closer 
to the magnetic pole and thus hotspot outer radius decreases. In 
case of a hypothetical circular spot, decreasing of the spot size 
p would result in increasing of the amplitude (expression 22 in 
IPoutanen & Gierliilskill2003 ). which contradicts the observations. 
A filled circular spot, however, is not what is seen in MHD simula- 
tions, instead a ring- or a crescent-shaped spot shapes are observed 
jRomanova et al 1 l2004h . In such a case, the emission is generated 
in a preferred sector situated closer to the disc. When the outer spot 
radius decreases, the preferred sector shifts closer to the magnetic 
pole of the star reducing the effective spot colatitude and leading to 
a decrease of the pulse amplitude as observed. 



5.3 Origin of the lag evolution 

The phase lags can be explained by a difference in emis- 
sivity patterns at differen t energies dOierlinski et al.l |2002| : 
IPoutanen & Gierliilskill2003h . At relatively high energies (above 
10 keV) Comptonization is dominating the emission. But as we 
proceed downwards from 10 to few keV, we observe gradually in- 
creasing contribution from the blackbody component which alters 
the emissivity pattern (towards more isotropic emission). There- 
fore maximum of flux is observed for different energies at dif- 
ferent pulse phases and this effect creates energy-dependent time 
lags. In agreement with this explanation, the observed lag saturates 
at the energy where blackbody component becomes insignificant. 
IGRI7511 demonstrates a smooth increase in the time lags of the 
fundamental during the outburst (Fig.|5) and a nearly constant time 
lags of the overtone (rather weakly constrained by observations). 
Several sources show that the lag does not remain constant during 
the outburst: in SAX JI808.4-3658 the lag value increases during 
the "sl ow decay", while in the "flaring tail" stage it starts to de- 
crease jHartmanet al]|2009h . In XTE J1814-338 there is also an 
obser vational evidence of the lag increasing in the end of the out- 
burst jWatts & StrohmaverllioOQ) . 

While the detailed modelling of actual pulse profiles is beyond 
the scope of this work, below we discuss the likely cause of the 
observed lag evolution. The pulse shape depends on many factors, 
such as the spot shape (which can be different for different emission 
components) and size, colatitud e of the magnet ic pole and the in- 
ner radius of accretion disc (see |Poutanenll2008l) . Changes in these 
parameters will affect the pulse shape, but not every of them can 
generate the observed increase of lags. 

To study the lag evolution, we employed the follow- 
ing toy model (fo llo wing th e framework develo pe d by 
Poutanen & Gierlinskil l2003l: IViironen & PoutanenI 12004 



rr 



-200 



Poutanen & Beloborodovl 20061 : llbragimov & PoutanenI 12009 
We have assumed a system inclination of 60° and a 1.4 M© 
star with radius of 10.3 km rotating at the pulsar frequency 
of 245 Hz. To model the emission anisotropy, we have cho- 
sen the angular dependence of emitting radiation in the form 
/(a) — Jo(l — hcosa), where h is the anisotropy parameter and 
a is the angle between the spot normal and photon direction in the 
spot frame. For the blackbody radiation h = 0, while for Comp- 
tonization in a slab of T homso n optical depth about unity h > 0.5 
I Viironen & PoutanenI |2004| : llbragimov. Zdziarski. & PoutanenI 
120071). and for the observed AMP pulse profi les h ~ 0.7-0 .8 is 
required (Poutanen & Gierlinski 2003; Pouta nen et al. I l2009l) . We 
generate two pulse profiles: one in the hard energy band, where 
only the Comptonization component contributes, and another at 
soft energies is a mixture of the blackbody and Comptonization 
components in equal proportions (we neglect changes of the 
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Figure 10. Time lag dependence on the anisotropy parameter h from 0.5 
to 0.7. The pulse maximum lags and the lags of the fundamental and the 
overtone ai'e shown as the solid, dotted and dashed curves, respectively. 
The dash-dotted line represents the pulse maximum lag between the black- 
body and the Comptonization component lightcurves. The contribution of 
two components to the soft X-ray band are assumed to be equal, the spot 
colatitude 9 was set to 15°, the powerlaw photon index of the continuum 
spectrum was assumed to be 1.9, and the angular radii of blackbody and 
Comptonization emitting spots were assumed 30° and 15°, respectively. 



Comptonization emi ssivity pattern w i th ene rgy, which is a good 
approximation, see llbragimov et alj |2007|) . We consider two 
geometries: a circular spot and a crescent spot, mimick ing the 
shape obtained in MHD modelling teomanova et alj2004h . 

For both assumed spot geometries, we have arrived to a sim- 
ilar conclusion. A rather small change of anisotropy parameter 
(from 0.5 to 0.6-0.7) turns out to be the best candidate for the lag 
evolution, as it can reproduce the observed lag increase of about 
~ 200/is (see Fig. llOt . On the other hand, a physically realistic 
change of the ratio between Comptonized and blackbody emission 
in the "soft" light curve, a change of the relative size of Comp- 
tonization and blackbody spot, a change in the spot size or a change 
of the effective spot colatitude were not able to produce lag evolu- 
tion large enough to match the observed value. 

Additionally, we can compare the aforementioned time lags 
with the lags between physical components (blackbody and Comp- 
tonization light curves). Very similarly to what is observed with the 
phase-resolved spectroscopy, we obtain the "component lags" that 
are nearly two times bigger than the lags computed as a mixture of 
two components. For the case shown on Fig. [Toj the "component 
lag" changes from —250 to —740 /is. 

Variation of h in principle could be related to a decreasing op- 
tical depth of the Comptonizing slab (see Table |2}, which would 
allow for more blackbody emission to contribute to the soft energy 
band. A decrease of tt from 1.85 to 1.74 would result in a '--^ 10 per 
cent larger blackbody flux, making it unlikely to account for the ob- 
served increase of the phase lag. To achieve a lag increase of 200 fis 
for a typical case shown in Fig.[TO| the ratio of the blackbody flux to 
the Comptonized one in the soft band should change by a factor of 
several, which contradicts the observations of a remarkably similar 
spectra seen during the outburst. This would mean that the intrin- 
sic angular emissivity of emitting plasma should not change much. 
However, the vertical structure of the accretion shock might change 
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with variations of the accretion rate causing changes in the emissiv- 
ity pattern. Thus, if the anisotropy parameter h is varying during the 
outburst, it might indicate variations in the accretion shock struc- 
ture. Alternatively, a physical displacement of the centroid of the 
Comptonization component (accretion shock) relative to the black- 
body emitting spot can cause the lag increase. 



6 CONCLUSIONS 

In this paper, we have analysed the data from IGR 175 11 obtained 
by RXTE and Swift during its September-October 2009 outburst. 
For the largest part of the outburst the source was in the "slow de- 
cay" outburst stage. However, after 24 days from the beginning of 
outburst, right before the source faded beyond RXTE detection limit 
(MJD 551 12), a very sharp decrease of the source flux was detected 
accompanied by a considerable phase shift of the pulse profile. This 
clearly indicates the end of the "slow decay" stage at that date. 

The energy spectra of IGR 175 11 can be well described with 
Comptonization in the accretion shock (Te ~ 30 keV, rx ~ 2) 
and soft ~1 keV blackbody emission originating from the NS sur- 
face. The SwiftfX.KT data allowed us to estimate the accretion disc 
parameters. We obtained the inner disc temperature of about 0.24 
keV and reasonable inner disc radius estimates located within the 
corotation radius. We also detected weak reflection and an iron line 
(EW ~ 60-100 eV). Spectral evolution was rather weak dur- 
ing the outburst. Only in the very end of the outburst (after MID 
55105), we detected a slight change in the optical depth tt of 
the Comptonized emission; such a spectral stability is known to 
be common for AMPs. 

The pulse profiles of the source were smooth and do not 
demonstrate any prominent narrow secondary features, which indi- 
cates that observed emission comes mainly from one emitting spot. 
The pulse harmonics demonstrate drop of their amplitudes during 
the course of outburst. The study of phase lag behaviour revealed 
a considerable increase of the lag from 150 to 400 us, an effect 
seen also in few other AMPs. A change in the anisotropy pattern 
of the Comptonized radiation could explain the observed time lag 
evolution, implying changes in the accretion column geometry or a 
physical displacement of the centroid of the accretion shock rela- 
tive to the blackbody spot. 
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